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A-1 PRIGIiANl I.TSTINGS 


1 . SCOPE 


This specitication establishes the dosiejn for a general- 
ized linear model analysis of variance (AOV) computer program. 

The requirements specification for this program was provided by 
the Research ^ Test, and Evaluation (RT&E) Branch of the Earth 
Observations Division (EOD) of The National Aeronautics and Space 
Administration/Lyndon B. Johnson Space Center (NASA/JSC) . 



/ 


2 . APPLICABLE DOCUMENTS 


The following documents ^ of exact issue shown r form a part of 
this specification to the extent specified herein: 

• "Generalized Linear Model AOV Program Requirements," 
provided by M. D. Pore/LEC. 

d Task Agreement 77-4, Job Order 81-127 


• TIRF 77-0042 


3 . SYSTEM DESCRIPTION 


3 . 1 HARDWARE DESCRIPTION 
Not applicable. 

3.2 SOFTWARE DESCRIPTION 

The purpose of the program GLMAOV is to implement an analysis of 
‘variance for experimental designs. The program is to complete an 
analysis of variance for unbalanced designs, designs with missing 
data and designs with multiple solutions. The program is in a 
general form and includes analyses of covariance and regression. 
Pseudo-inverses of the model matrices required in the analyses 
will be computed using an algorithm identified as "Greville's 
Method . •' 

The program is coded in the Univac 1108/EXEC 2 Fortran V language 
but using only Fortran IV-G capability for future implementation 
on the IBM/370. 

3.2.1 SOFTWARE COMPONENT NO. 1 (GLMAOV) 

GLMAOV functions as a driver program for the generalized Linear 
Model Analysis of Variance Program (LMAOV) . It reads in the 
number of analyses to be taken, the dimensions of the input math 
models for each analysis and then computes base addresses for 
arrays used in LMAOV. The general linear model used is: 

Y = X3 + e 

with 3 subjected to the constraint: 

R3 = t 


The hypothesis tested is 

H^: X3 = h 



against the alternative 

X0 ^ h 


3. 2. 1.1 Linkages 

GLMAOV is not referenced by any other program. It calls sub 
routine lmaov. 


3.2. 1 . 2 Interfaces 

Interface between GLMAOV and LMAOV is done via the calling argu- 
ments of LMAOV. 


3 . 2 . 1 . 3 Inputs 
GLMAOV requires an 

Column 

1st card 1-3 
5-6 
8-9 
11-12 
14-15 
17-26 


input card deck as follows: 

[Right justify all numbers] 

Description 
Number of analyses to be taken 
NP (Row dimension of Y and X) 

NQ (Column dimension of X, R, and X) 
NM (Row dimension of R and t) 

NN (Row dimension of X and h) 

TOLENC (Tolerance Level) 


Subsequent 

Cards 


Required inputs for first analysis for subroutine 
LMAOV; see Software Component 2 (LMAOV) , section 

3.2.2 


1-2 NEW NN (Row dimension of next set of X and L) 
[inputs for second analysis for subroutine LMAOX^ 




Column 


Description 

[ Inputs for third analysis 
for subroutine LMAOV 


End-of-File card 

3. 2. 1. 4 Outputs 

If the math models used are found to be too large dimensionally 
for* the space allocated, GLMAOV outputs this message and termi- 
nates : 

****** Dimensions of math models are too large for use in this 
program ** Storage capacity of GLMAOV is exceeded. 

3. 2. 1.5 Storage Requirements 
To *be determined. 

3 . 2 . 1 . 6 Description 

GLMAOV first reads in the number of analyses to be taken, the 
dimensions of the math models [See 3.2.2 SOFTWARE COMPONENT 2 
(LMAOV)] and an input tolerance level. If TOLENC is not input, 
TOLENC defaults to ,0001. GLMAOV then computes base addresses 
of all arrays used in LMAOV from the dimensions that were read. 
GLMAOV has a large working array of size 42,000 which will be 
divided using the base addresses that were computed. These base 
addresses are passed to LMAOV via the calling arguments as start 
ing addresses of working arrays for LMAOV. The call for LMAOV 
is in a DO- loop from 1 to number of analyses. A dimension for NN 
must be input before each analysis whether NN will be changed 





or not 




3. 2. 1.7 Flow Chart 


Not applicable. 

3 . 2 , 1 . 8 Listing 
Sec Appendix A 

3.2.2 SOFTWARE COMPONENT NO. 2 (LMAOV) 

LMAOV is a generalized Linear Model Analysis of Variance Program 
for experimental designs. It will complete the analysis for un- 
balanced designs > designs with missing data and designs with 
multiple solutions. The model is written in a general fo.rm and 
includes analyses of covariance and regression. The general 
modal used is: 

Y = xa + e 

with a subjected to the constraint: 

R$ = t 


The hypothesis tested is 

xa = h 

against the alternative 

xa h 

3.2c 2 . 1 Linkages 

LMAOV is called from GLMAOV and will reference TRANSP, MULTMX, 
SUBMX and PSINV. 

3 . 2 . 2 . 2 Interfaces 

Interfaces provided by calling sequences. 




3. 2. 2. 3 Input! 

Calling Sequence; 

Call LMAOV (YVEC, XMX# BVEC, RMX# CTVEC, RINV, TVEC, GMX/ GINV, 


HINV, CHV, HVEC, WKR, 
WKRS, CIRV, MMX, ZWRK, 
WORXST, NP, NQ, NM, NN 

XTRP, HMX, GWKR, CHVEC, RWKR, CMX, CHMX, 
ZWST, PWRK, MINV, PWST, ZTRP, PPWRK, EWRK, 
, NALYS, TOIENC, NWS) 

a 

Parameter 

Dimensions 

In/Out 

Description 

YVEC 

NP 

In 

Random vector of observations 
(Y vector) 

XMX 

(NP, NQ) 

In 

Known matrix that includes the 
design and regression matrices 
(X' matrix) 

BVEC 

NQ 

In 

Vector of random and non-random 


parameters including higher 
power terms and interaction 
terms ($ vector). The test 
variables, regression varia- 
bles and con-commitant varia- 
bles are subsets of 0 


RMX 

(NM, NQ) 

In 

R matrix 

CTVEC 

(NM) 

In 

Working array 

RINV 

(NQ, NM) 

In 

Working array 

TVEC 

(NM) 

In 

t vector 

GMX 

(NN, NQ) 

In 

\ matrix 

GINV 

(NQ, NN) 

In 

Working array 

HINV 

(NQ, NN) 

In 

Working array 

CHV 

(NQ) 

In 

Working array 

HVEC 

(NN) 

In 

h vector 
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p«raiMt«r 

DiitMinaion 

In/Out 

Description 

WXR 

(NO, NQ) 

In 

'^Miking array 

XTRP 

(NQ, NP) 

In 

Working array 

HMX 

(NN, NQ) 

In 

Working array 

GWKR 

(NN, NN) 

In 

Working array 

CHVEC 

(NN) 

li. 

Working array 

RWKR 

(NM, NM) 

In 

Working array 

CMX 

(NQ, NQ) 

In 

Working array 

CHMX 

(NN, NM) 

In 

Working array 

WKRS 

(NQ, NQ) 

In 

Working array 

CINV 

(NQ, NQ) 

In 

Working array 

MMX 

(NP, NQ) 

In 

Working array 

ZWRK 

(NP) 

In 

working array 

ZWST 

(NP) 

In 

Working array 

PWRK 

(NQ, NP) 

In 

Working array 

MINV 

(NQ, NP) 

In 

Working array 

PWST 

(NP, NP) 

In 

Working array 

ZTRP 

(NP) 

In 

Working array 

PPWRK 

(NP, MP) 

In 

Working array 

EWRK 

(NQ, NQ) 

In 

Working array 

WORKST 

1 

In 

Starting address of the remain- 
ing storage of the large working 
array in GLMAOV 

NP 

1 

In 

Row dimension of Y and X 

NQ 

1 

In 

Column dimension of X, R, and 
•X 




Parameter 

Dimension 

In/Out 

Description 

NM 

1 

In 

Row dimension of R and t 

NN 


In 

Row dimension of X and h 

NALYS 

1 

In 

Number of analyses to be taken 

TOLENC 

1 

In 

User input tolerance level 

ms 

1 

In 

Last computed base address 
from GLMAOV 


Xnput of the elements of the known vectors and matrices are read 
from cards in the following manner i 

1. X(XMX) and y(YVEC) are read simultaneously. First one row 
of X is read followed by the first element of Y. The second 
row of X is read followed by the second element of Y, etc. 
Hence NP lines of NQ + 1 elements. 

2. R(KMX) and t(TVEC) are read in the same manner as (l.)y 
therefore NM lines of HQ + 1 elements. 

3. X(GMX) and h{HVEC) are also read in the same manner as (1.) 
and (2,)y hence NN lines of NQ + 1 terms. 

For a row of terms, a maximum of 8 terms may be input per card, 
the first term in columns 1-10, second in columns 11-20, third 
in columns 21-30, etc. until (as for X) NQ terms are read from 
as many cards as needed for NQ terms. The card following will 
contain the first element of Y (columns 1-10) then the same 
process is repeated for the next row of X followed by the next 
term of Y and so on. Cards for R and t follow in the same 
manner; then cards for X and h. If more than one analysis is 
parfornad, the data for R and t and the data for X and h would 
follow in the saaie fashion as before. In the case that more than 
one analysis is taken, a card (punched i;n col, lr2). fs necessary 
that contains the row dimension (NN) of X and L. This is done so 
that the number of rows of X and L may vary from analysis to 
analysis. Even if the dimension remains the same with sucessive 
analysis, NN must be input with each set of inputs. 

f 


The user must use care in having the correction number of terms 
for each row for the correct number of rows. An example of the deck 
setup is located on the next page. 



DECK SETUP FOR PROGRAM, 




HH 


Decks of 
R and t 
■Smd X and 


G.GG GGG.G GG. 


GG.GG 


H.H 


h follow 
for N nuir.‘ 
ber of 
analysis 


GG.GG GG.G GGG.G GG.GG 


N analyses 
2nd analysis 




T.TTT 


c 




RR RRR.R RR.R 














m 


New row dimension 
of X 


HH 


1st term of h 


GGG.GG GGG.G GG.G GGGG.G 1st pw of X 

Last term of t 


TT . T . TTT 


RR.RR RR.R RRR.R RRRR.R 2nd row of R 

1st term of T 


TT.T 




r 


R.RRR R.RRR RRR.R RRRR. R 1st row of R 

W'TTT 


Last term of y 


XXXX.XX XXXX.XX XXXX.X XXX. X 2nd row of X 

1st term of y 


t analysis / ^ 

I ^ • X 

^^XX.XX xkxxx.x XXXY. xxxx. xxx.xxx 






C.l 


n 


21 


31 41 


80 1st row of X 


NALYS 


NP 


C.l 


NO 

3 


NM NN TOLENC 

n 14 17 


(Front of Deck) 


Also, input terms are presumed to be decimal (floating point) 
numbers of a maximum of 8 digits with a maximum of 6 digits to 
the right of the decimal and should be located within a card 
field of 10 columns, starting at column 1, i.e. , 1-10, 11-20, 

21-30, 31-40, 71-80. 

3 . 2 . 2 , 4 Outputs 

The following messages are printed if computations described in 
Section 3. 2, 2. 6 are found to be true* 

2^^ if t it it it ie if if if 

REJECT Hq with PROB (TYPE I ERROR) = PROB (TYPE II ERROR) = 0. 
THE RESTRICTIONS AND HYPOTHESES ARE SEPARATELY CONSISTENT, 

BUT NOT JOINTLY COMPATIBLE. 

2 , *****#*»»** 

Hq is a legitimate hypothesis BUT THE MODEL IS NOT MATHEMAT- 
ICALLY CONSISTENT. 

RECONSIDER THE ASSUMED RESTRICTIONS. 

3 ^ ********** 

REJECT H^, WITH PROB (TYPE I ERROR) = PROB (TYPE II ERROR) = 0. 
THE ASSUMED RESTRICTIONS ARE MATEHMATICALLY CONSISTENT, BUT 
THE HYPOTHESIS IS NOT CONSISTENT. 

4, ********** 

REJECT Hq with PROB (TYPE I ERROR) - PROB (TYPE II ERROR) = 0. 
NEITHER THE ASSUMED RESTRICTIONS NOR THE HYPOTHESIS IS MATHE- 
MATICALLY CONSISTENT. 

5 , ********** 

THE HYPOTHESIS FUNCTION, 3 IS NOT ESTIMABLE. ANALYSIS CON- 
TINUES. 

If 1. through 4. of the above output is not printed, LMAOV prints 
an Analysis of Variance (AOV) table, the Best Linear Estimate (BLE) 
unrestricted with the covariance matrix and squared multiple 


correlation coefficient, the BLE restricted only by the model 
restrictions vrith the covariance matrix, and the BLE restricted 
by the hyposthesis with covariance matrix and squared multiple 
correlation coefficient. Refer to section 3. 2. 2. 6 for definition 
of values of the AOV Table, covariance matrices, and squared 
multiple correlation coefficients. 


AOV TABLE 


Source of 

Degree of 

Sum of 

Mean 


Value of 

variation 

freedom 

squares 

squares 

P • 

- statistic 

Due to 3 

^1 

Q3 




Due to 3« 
(Unad j . ) ” 


Q.. 





4 



NP 


Due to 3jj 

R 1 -R 2 

Qi 

MSI 

p = gi 

(ADJ) 





1 2 

Error 

P - R^ 

Qo 

MS 2 



Total 

P 

0 





The Best Linear Estimate (BLK) unrestricted is; 

[r] 

with covariance matrix; 

[GOV (f)] 

And squared multiple correlation coefficient; 

r2 

The BLE restricted only by the model restricted is; 

[a] 

with covariance matrix 

Ccov S] 

The BLE restricted by the hypothesis is; 

ci] 

with covariance matrix 

[COV §] 






And squared multiple correlation coefficient: 


LMAOV also outputs an error message and 
storage is exceeeded: 

* * * * * dimensions of math models are 

GLMAOV * * STORAGE CAPACITY OF LMAOV IS 

3 . 2 . 2 . 5 STORAGE REQUIREMENTS 
To be determined. 

3. 2. 2. 6 Description 

Initially LMAOV reads input data. Y(YVEC) and X(XMX) are first 
read simultaneously in a nested DO- loop. Each row of input is 
a row of X followed by the respective element of Y so that there 
are NP rows of NQ 1 terms. The same operation is done for 
R(RMX) and t(TVEC), hence NM rows of NQ 1 terms, and X(GMX) 
and h(HVEC) , NN rows of NQ + 1 terms. 

The computations in LMAOV require the subroutine PSINV which 
computes the Moore-Penrose pseudo- inverse of a matrix. Given 
any non-zero ^ x qj matrix A, the pseudo- inverse is the unique 
matrix A'*’. (Refer to section 3.2.3 SOFTWARE COMPONENT NO. 3 
(PSINV)). 


terminates execution when 

TOO LARGE FOR USE IN 
EXCEEDED * * * 


H » X(I - R"^R) 



Initial computation: 


I 

'fT 


A t«st is made to determine if h(HVEC) and t(TVEC) are zero 
actors. Comparison to zero is actually a test against an arbi- 
trarily small value, greater than absolute zero. The test that 
is made is: 

2 (a^**2) < (T0LENC**2) * AMAX (l , E (,a^**2)) 

Zf true, the vector is considered to be a zero vector (if TOLENC 
is not input by the user, default TOLENC is .0001). 

Zf both h and t are zero vectors, LMAOV computes values for the 
output AOV table. Zf either/both h and t are non-zero vectors, 
the following tests are made: 

A test for XX^h ■ h and RR'*’t = t is made. Equality of vectors 
is determined by: 

■ ABS(Z((a^ - b^)**2)) < (TOLENC**2) * AMAX |l , (e (a^**2) + E(b^**2))j 

Zf the test is true, the vectors are considered equal. 

1. If (XX^h ■ h> and (RR'^t ^ t) the following message is written 

and execution cn this set of data is terminated: 

* * * * 

Hq is a legitimate hypothesis but THE MODEL IS NOT MATHEMATI- 
CALLY CONSISTENT. RECONSIDER THE ASSUMED RESTRICTIONS. 

4 * 4 * 

2. If (XX h h) and (RR t - t) the following message is written 

and execution on this set of data is terminated: 

* * * * 

REJECT Hq with PROB (TYPE I ERROR) = PROB (TYPE II ERROR) = 0. 
THE ASSUMED RESTRICTIONS ARE MATHEMATICALLY CONSISTENT, BUT 
THE HYPOTHESIS IS NOT CONSISTENT. 




+ + 

3. If (XX h h) and (RR t ^ t) the following message is written 
and execution on this set of data is terminated: 

****** ' 

REJECT Hq with PROB (TYPE I ERROR) = PROB (TYPE II ERROR) 0 
NEITHER THE ASSUMED RESTRICTIONS NOR THE HYPOTHESIS IS 
MATHEMATICALLY CONSISTENT. 

4. If ( XX^h a h) and (RR"*’t = t) another test is made: 

HH"*‘(h - XR'^t) - (h - XR'^t) 

If the above test is true, LMAOV conducts the test in 5. and 
computes the values for the AOV table; if not true, the 
following message is written and execution on this set of 
data is terminated; 

****** 

REJECT Hq with PROB (TYPE I ERROR) = PROB (TYPE II ERROR) = 0 
THE RESTRICTIONS AND HYPOTHESES ARE SEPARATELY CONSISTENT, 

BUT NOT JOINTLY COMPATABLE. 

5. If H[X(I-r'''r) ] [X(I“R’^R) H 

then print: 

THE HYPOTHESIS FUNCTION, IS NOT ESTIMABLE. ANALYSIS 

CONTINUES . 




COMPUTATIONS FOR AOV TABLE AND BEST LINEAR ESTIMATES 
Calculate : 

M - X(I - R'^R - 
C » (I - R'^R)X^X(I - r'*‘r) 

Z - Y - XR'^'t - XH‘^(h - XR'^t) 

P - length of Y » NP 

B p - Tr(I - XC X ) rounded and truncated to an integer 

Rj “ “ Tr (XC‘*’x^’'>- MM^) rounded and truncated to an integer 

T + T 

Qg « Z^XC X Z 
T + 

Qj “ Z MM Z 

Z^(XcV - MM*^)Z 

Qq = Z^(I - XC’^X^)2 
Q =* Z^Z 



The following AOV tcible is printed from the above computations. 



Source of 
variation 

Degree of 
freedom 

Sum of 
squares 

Mean 

squares 

Value of 
P - statistic 

Due to B 

«1 




Due to B(N) 
(Unad j * ) 

R 2 

Q2 



Due to 0(H) 
(Adj . ) 

«1 " ^2 

Qi 

MSI 

F 

Error 

- 

Qo 

MS 2 


Total 

P 

Q 




The the following values are calculated: 

YMEAN » eY(i) , i * 1 thru NP 
NP 

+ 

r = X Y 

cov(r) s (x’^x)'^ 

^ ^ ” r X Y ^ NP (YMEAN) 

y'^Y - NP (YMEAN) ^ 
a = R t + C^X-^ (Y - XR t) 

COV(a) = 

J = R'*'t + H‘^(h - XR'^'t) + M"^[Y - XR'*'t - XH'^(h - AR'^'t) ] 

COV(§) = (m'^m)'*' 

= [ 7 !^ (MB ) - NP ( YMEAN^ ) ] ^ 

( Q - NP ( YMEAN^ ) ) (0 '^M^M P -NP ( YMEAN^ ) ) 

The above calculations are printed as such: 

The Best Linear Estimate (ELSE) unrestricted is: 

[r] 

with covariance matrix: 

COV ( f ) 


and squared multiple correlation coefficient: 

[rR=l 


Th« B«st Linear Estimate restricted by the model restrictions is: 

[3l 

with c^vVariance matrix: 

COV (o) 

The Best Linear Estimate restricted by the hypothesis is: 

m 

with covariance matrix: 

COV (g) 

and squared multiple correlation coefficient: 

m^] 


All transposition operations are done in subroutine TRANSP, matrix 
multiplication in subroutine MULTMX, and matrix subtraction in 
subroutine SUBMX. 

All storage and working storage for LMAOV is passed by address 
into LMAOV via the calling argument from 6LMA0V. This was done 
to provide the user with the maximum possible storage available 
for batch runs on the UNIVAC 1108. If the user exceeds this 
allocation of space / a message will be printed stating the vio- 
lation. In which case the user should readjust the math model. 

LMAOV uses the variable WORKST as the starting address of the 
remaining working space from ARRAY in GLMAOV. LMAOV will compute 
all base addresses for arrays in PSINV and divide WORKST appro- 
priately. If storage space is exceeded, a message is written 
and the program stops execution on that particular set of data. 

jpK 


3 . 2 . 2 . 7 Flow Chart 


Mot applicable. 

3 . 2 , 2 . 8 Listing 
F^e V^pendix A 

3.2,3 SOFTWARE COMPONENT MO, 3 (PSINV) 

The function of PSINV is to compute the Moore-Penrose pseudo 
inverse of a matrix. The defining properties of the pseudo- 
inverse are as follows: Given any non-zero ^xg] matrix A, the 

pseudo-inverse is the unique matrix^ denoted A*^ , such that 

1. AA‘‘"a * A 

2. A'^’AA'*' « A**" 

3. (AA’*')'^ » AA'*’ 

4. (A‘‘‘a)'^ » A'^'A 

If A is a square and full rank matrix, then A has an inverse and 
a'*’ « a“^. If a is a non-zero|]pxl] vector then 

a'*’ » (a'^a)"^a'^. 

IT 

(A denotes the transpose of the matrix A. ) The algorithm 
used for this computation is identified as "Greville's Method." 

3. 2.3.1 Linkages 

Subroutine PSINV is call«?d from LMAOV. PSINV calls TRANSP, 
MULTMX . and SUBMX. 

3. 2. 3. 2 Interfaces 

Interface between PSINV and all other subroutines is via 
the calling arguments. 


3pi^ 


3 . 2 . 3 . 3 Inputs 
Calling sequence is: 

Call PSINV (AMX, OUTMX, NROW, NCOL, TOLENC, FCOL, TOOL, DK, CK, 
AK, DTR, BK, SCRACH) 


Parameter 

Dimension 

In/Out 

Description 

AMX 

(NROW, NCOL)’ 

In 

Real matrix from which 
pseudo-inverse will be taken 

OUTMX 

(NCOL, NROW) 

Out 

Real matrix which will be the 
pseudo- inverse matrix of AMX 

NROW 

1 

In 

Row dimension of AMX 

NCOL 

1 

In 

Column dimension of AMX 

TOLENC 

1 

X 

In 

Tolerance level 

FCOL 

('^ROW) 

In 

Working array 

TCOL 

(NROW) 

In 

Working array 

DK 

(NCOL) 

In 

Working array 

CK 

(NROW) 

In 

Working array 

AK 

(NROW) 

In 

Working array 

DTR 

(NCOL) 

In 

Working array 

BK 

(NROW) 

In 

Working array 

SCRACH 

(NCOL, NROW) 

In 

Working array 


3 . 2 . 3 . 4 Outputs 

PSINV outputs OUTMX matrix. OUTMX is the pseudo-inverse of AMX 
matrix. 




3. 2*3. 5 Storage R<quir«menta 
TO be determined. 


3 . 2 . 3 • 6 Description 

PSINV is a subroutine that uses *Greville*s method' as the algo- 

4 * 

rithm to compute the pseudo- inverse of a matrix, denoted A . The 
matrix to be inverted and its dimensions are input parameters. 

Computing the Pseudo-Inverse 

Computation is accomplished by taking one column of the input 
matrix ^ time to produce one row of its inverse [bUTMxJ . 

a^ will be the first column of AMX. 

If e(a^**2) < (T0LENC**2) * AMAX(1, e (a^ **2) 

then a^ is considered equal to zero. (Test for zero vector) . 

If » 0 then A^ ■ 0 (row vector, first row of OUTMX) . 

^TOLBNC is input from LMAOV as a testing level for determining 
the closeness of the vector to a zero vector] . 

+ + T “"I T 

If aj^ / 0, then 

T ••X 

Note that (a^^a^^) results in matrix multiplication of (1/NN)X 
(NN, 1) so that a (IXI) SCALAR results. Consequently the inverse 
is taken as 1. /SCALAR. 

AMX is then partitioned by columns in a DO-loop from K = 2, NCOL. 
A„ is the sub-matrix of AMX consisting of K columns and a„ is the 

K A 

column. 





N = K - 1 


thfn d|^ is coii|»ut«d «• 




then ie computed as 


=K “ »K - Vk 


then is determined as 

N - =!i' =K “ 


or 


^0? + 

N 


'K 


1 + 


m I 

dj^dn 


if C 


K 


Resulting In as 


b„ being an added row after each Iteration until K » NCOL, 

A 

All transpositions are done by subroutine TRANSP, matrix multi- 
plication by subroutine MULTMX, ^nd matrix subtraction by sub- 
routine SUBMX. 

All working array storage is taken from already existing storage 
in LMAOV where base addresses are computed accordinij to the 
dimensions of AMX and OUTMX. 






3.2.3. 7 Plow Chart 


Not applicable. 

¥ 

3 . 2 . 3 . 8 Listing 
See i^pendix A 

3.2.4 SOP^t^ARE COMPONENT NO. 4 (TRANSP) 

The function of TRANSP is to transpose a matrix of real elements. 

3.2. 4.1 Linkages 

TRANSP is referenced by LMAOV and PSINV. TRANSP does not require 
any other subroutines. 

3 . 2 . 4 . 2 Interfaces 

Interface between TRANSP and the two programs that call it are 

via the calling arguments of TRANSP. 

¥ 

3 . 2 . 4 . 3 Inputs 
Calling Sequence: 


Call TRANSP 

(A, B, M, N, 

MX, NX) 


Parameter 

Dimension 

In/Out 

Definition 

A 

(M, N) 

In 

The matrix to be transposed 

B 

(N, M) 

Out 

The transpose of A matrix 

M 

1 

In 

Row dimension of A 

N 

1 

In 

Column dimension of A 

MX 

1 

In 

First dimension of A as 
specified in the DIMENSION 
statement of the calling 
program 



Parameter 


Dimension 


In/Out 


Definition 


NX 1 In First dimension of B as 

specified in the DIMENSION 
statement of the calling 
program 

3. 2.4.4 Outputs 

B matrix is output as the transpose of A. 

3 . 2 . 4. 5 Storage Requirements 
To be determined. 

3. 2. 4.6 Description 

The transpose of A (M, 
elements are: 

3. 2. 4. 7 Flow Chart 
Not applicable. 

3.2. 4.8 Listing 
See y^pendix A 

3.2.5 SOFTWARE COMPONENT NO. 5 (MULTMX) 

MULTMX multiplies two matrices of type REAL. 

3. 2.5.1 Linkages 

MULTMX is referenced by LMAOV and PSINV. MULTMX does not require 
any other subroutine. 


N) matrix is B (N, M) matrix whose 


B . . =s A. . 
31 13 




3 , 2 . 5 . 2 Interfaces 


Interface between MULTMX and the two subroutines that call it is 
via the calling arguments of MULTMX. 


3.2.5. 3 Inputs 
Calling Sequence: 

Call MULTMX (A, B, C, M, N, K, MX, NX) 


Parameter Dimension 
A (M, N) 


B (N, K) 


C (M, K) 


In/Out Definition 

In Two dimensional array 

containing elements of 
multiplicand matrix 

In Two dimensional array 

containing elements of 
multiplier matrix 

Out Two dimensional array 

containing product of A 
and B 


M 1 

N 1 

K 1 

MX 1 


In First dimension of A and C 

In Second dimension of A and 

first dimension of B 


In Second dimension of B and 

C 


In First dimension of A as 

specified by DIMENSION 
statement in the calling 
program 


In 


^7 


First dimension of B as 
specified by DIMENSION 
statement in calling 
program 


NX 


1 


3 . 2 . 5 . 4 Outputs 

C matrix (M# K) which will be the product of A matrix (M, N) and 
B matrix (N, K) . 

3. 2. 5. 5 Storage Requirements 
To be determined. 

3. 2. 5. 6 Description 

The product of the (M, N) matrix A and (N, K) matrix B is a (M, K) 
matrix C whose elements are defined as: 

N 

C . . = \ '' a . X b . 

13 2-1 IL hi 

L=1 


3. 2. 5. 7 Flow Chart 
Not applicable. 

3. 2. 5. 8 Listing 
See Appendix A 

3.2.6 SOFTWARE COMPONENT NO. 6 (SUBMX) 

SUBMX subtracts two matrices of type REAL. 

3 . 2 . 6 . 1 Linkages 

SUBMX is referenced by LMAOV and PSINV. SUBMX does not call any 
other subprograms. 

3. 2. 6. 2 Interfaces 


Interface between SUBMX and LMAOV and PSINV is only via the 
calling arguments of SUBMX. 


3. 2. 6. 3 Inpttti 

Calling Sequaneei 

CALL SUBMX (A, B, C, H, N, MX) 


Parameter 

Dimension 

In/Out 

Definition 

A 

(M, M) 

In 

Two dimensional subtrahend 
matrix 

B 

(Mr N) 

In 

Two dimensional minuend 
matrix 

C 

(Mr M) 

Out 

Two dimensional matrix 
containing difference of 
A and B 

M 

1 

In 

Row dimension of Ar B, and 
C 

N 

1 

In 

Column dimension of Ar B, 
and C 

MX 

1 

In 

First dimension of Ar B, 
and C as specified in 
DIMENSION statement of the 
calling program 


3. 2. 6.4 Outputs 

Matrix C (M, N) will contain the difference of matrix A (M, N) 
and matrix B (M, N) . 

3. 2. 6. 5 Storage Requirements 
To be determined. 

3. 2. 6. 6 Description 

SUBMX consists of two nested Do- loops where I = 1, N and J = 1, 
M. The difference of the elements of the two matrices where 



C ■ A " and B " b^j iat 

c . j ■ a . . - b . . 

a.j ij 13 

3. 2. 6. 7 Flow Chart 
Not applicable. 

3. 2. 6. 8 Listing 
So6 /^pendix A 



4 . OPERATION 


4.1 OSER DOCUMENTATION 
Run Deck Set up: 

@ Run Badge ID, Div CODE, BOX #., Proj#, Proj#, C, Time, Page 

@ SCH 7T = 1 

@ ASG A » X01189 or X01180 

@ XQT CUR 

TRW A 
IN A 

@ XQT GLMAOV 

[ Input data as 

described in Section 3.2,1, 3 and Section 3. 2. 2. 3 
@ FIN 




Col 61 
NAME 


APPENDIX A 


PROGI^W LISTINGS 






9 FOR SOBhXiSUBHK 

tSls'UJSL.'Sgras'JSlSH'jfe'S'.l; 

subroutine SOBHX tNiHy point nc*376 

storage used: CODEII} CCOIII; OATAlr-i sceos^j BtAKiC C0NNCNI2, C3O0C0 
^ external reference^ iBuOCfc, NAhEI 

D003 NERK3S 

storage aBSISNHENT IBLOCK, TTPE, RvlATU-E LOCATlCNt t»A«El 
0001 000037 loSs QOOI ifCcRC Use 0003 1 CLECql I 


C^lOO 

!• 

C* 

qmcc 

2^ ' 

C« 

Odioo 

3* 

c* 

OCiOi 



OOIO> 

S» 


Coio^ 

6* 


05107 

7» 


Cri 1 z 



CC1I3 

9» 

to 

Oci 16 

ID« 


00117 

!!• 



SbbKOUTlNF SUBmX I’OtS •*£*(. NATHIX SLbTKA^TlON 

SUBROUTINE 5UBhX‘< ihtt»»intNX> 

OiNENSIUN AIhX.nIi Bln^tNlt LIRX.Ni 
Ut 10 JBI.N 
OO 10 1 * 1 . H 
C(1 •J)>All .Jl- I • J) 

Continue 

heturu 

End 


CHO or COHFlLATlnh: 


NO i.I^CNoST^CS, 




25-HAR tB 




8 FOR HoUTHAiHUI.TKX 

^O'^TRAH V E»Ec II UEVrt 2i.ii ,-«E*EVfc LEtftL EIZOIDBIO^I 

This cohpiuation »*s uohE oh »s m*r 7b *t 10:3^:29 

SUOROUTlNt HOWTHA EwTHt PoImT 0C»‘132 • 

storage USEOj COQElII OCOlSGi UATAtO BUAHK COMHUNIZJ GqoOCO 

EXTERMAI. REFERENCES IbuOCX, HAhEJ 
UDU3 NERK3S 

storage aSSICNHENT cBLOcX, TYPt. KfEATIvE GoC^TIOMi NaHE> 


0001 OOOOSQ 10*6 
0300 I 030002 J 


6001 290r&S Hit 

0000 I ootrni** K 


ooai oooos? 11S& 

003u 0 033GUU 4 


0300 1 000003 1 


Oosa BOOQO* 1N3PS 


OnlQO 

i« 

C* 



Doloa 

2 » 

C* 

MUUThX SUhKUUTinE Dots HAtKiX HULT 1 PLVLaTi ON IN DOiiBt£ PHtCiSlUN 


OOIOQ 

3* 

C* 

qOT RESToHFS KtSUUT 0^ ^lUiiLE PR£Cl2^I0w TROHC^Ti-Or 


OOI QO 

H* 

C* 



OQiav 

b» 


iUBHOUTiNE 


00103 



UlHEflSlON X4HX.L** B ^ i N ) # “ C < «X » H 1 


coios 

/• 


OUUBLE PRrLlSI*)H * 


ODIOS 

80 


'Do in hN 


ODUO 

90 


Do ID i«l»N 


CO 1 1 3 

IV 


Cl I 


00 i 1 H 

1 1* 


DO 10 K»I.L 


00117 

12 * 


UIA|^> 


OoUO 

’ 13* 


Cil»J) « to * 


00121 

IS* 

10 

CoNT i f|UE 


OoUi* 

IS* 


KCTURN 


Co 126 

46« 


ChD 


Si 

tND OF 

COUP I EAT Ii»H; HO rj I AG'«nSTi CS. 


V 

V 



- 



t 


i 


M fOH TRaNSP, TRAH5P 

FUKTRaU V ExEc II LEVfL 75a 
This COHPiUATlON iSTAS done on >5 HAK 7b aT 


iCAE^a tEvEt eizoiQUio^i 


I 

25 HaR 7 * 10331 { 30 * 3*2 ’ 


suartouTiNt transp emTry point rjonaT^j 

storage u^EDS codec I) 11001 o6i OA|'A(d CD0»;3^r rUaNK COmKOUT2| OooQoo “ * . . I 


EXTERHaL REFEREHCES fbLOCK* NAHt) 

0003 NERN35 ^ 

«r 

STORAGE ASSlGNHLift iBuoCKt TTpt. RpL^TIvE 1--oC*TIon* NAHT > " - 

oDol oooo'sz loSc ono» siaoH-^ Jioo Bns3 i oosaoo i ooos ooooo> lajpt oooo i oooooi j 

* ‘ 

iloo I* c» : 


J 


□ OlOO 
SOlOO 
QDlOl 
00103 
OdIOR 
C01O7 
OoU2 
0 q |13 
Do I 16 
00117 


2 « 

3« 

Hm 

Sm 

6 « 

7* 

8 * 

9 ^ 

ID* 

II* 


C* 

c* 


IQ 


Transp does real hatki^ transposition 

SUBhOUTiUp TRANSFl \ , B . *» , » ,hx iNA I 
0 1 HENS I ON AChx*N^ * 

00 IP 

Do 10 

be J»I l«At] *J| 

Continue 

return 

End 


END OF COKPlUATl oN: NO ** I AOUoST ^ CS. 



1 




1 

. J 


- ■ 


iN rrx 


8 FOR fSInV.PSINV , 25"|iA'* »• 

UHlVfcC m»6 FOKTRAN V Cttc II LEVrU 25i ,*«EAtts tEvEL El20lC0l0*l 

This cohpiuation kas done on je n*h 7b »i 


SuBKoutINE PSJNy EwlRr POl*.r . 

storage used; cooiiii 00057*1 ; oataidi blank cokhonizi ooococ 


ekternau references IOeOCK. hanei 

0D03 TRANSP 
ODOR HUUMX 

.0005 SOBHA 

0CU6 NERK3S 


SToRACE ASSUNKEnI IBeOCK. TTpE, RkLATIvE ‘•oCaTIOn* NANE* 


OCOI OGDOlS 10*6 
OOQi OCOlSI 1526 
OCOI ODOOIT 3CU 
0000 R 300007 ON 
0000 1 ODOEfOS N 


OCOI 
0001 
0001 
SOQO I 

0030 R 


Jic-.H** 117* 

:i^Z72.\ I6Ht, 
133 I 2f *3^ 

iSEr.'X Scalar 


0001 

cool 

cosi 

5CS3 
CDS 3 


CLClSj 12*6 
CCOJC? 20lfc 
tyOJlJ yoL 
rtcoij inJp* 
OCCCCC aA 


gggiU IHI SIS! 18:I?$ nil 

QBQJSZ VSU foaa R SB300Z CC 

I coooo^ ic doss I osooo* 


0*8 100 

act 00 

OSIOD 

OOlOD 

CCIOC 

Do too 

SOlOO 
Co 100 
OClOO 

(01 Solo? 

^ OOIDI 
iK C0I03 
C0103 
00103 
00103 
OOlO«f 
00105 
001 10 
^ Dot II 
00113 
. COll^ 
ecii^ 
001 
C012I 
00122 
0CI2H 
0012H 


I« 

2R 

3» 

?• 

OR 


I?: 

I!: 

|H* 

15* 

17* 

tae 

IV* 

20* 

21 * 

22 * 

23* 

ZH9 

25* 

2 ** 

27* 

2a* 


C* 

C* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 


c* 

c* 


10 

c* 

20 

c* 


Subroutine psikv cokpu^es thE HoOHE*PEfcR2^t Pstuoo-iNVtKSt of a 

NATRlXf otNOTtri A^l 05*N<» AtOoHllMN aS tKtViUuE^S NLThOO* 

If *A* SoUA»'E *K0 full KAHK T«EN •A* AN 1NVLrs£ 

also If •a**is*^ nonI^ho pu vectuk tke*^ 

e« « I 


B«tI Be1h6 transfuse of a 

ANO t-II UfNut|F-<« iNVfcKSt, 


SUBROUUNr PsIa.vU«*# 0^TN» •PH0N,NC0ttl0tFNCffCUL»tC0LiCK»CR*AR» 

*UlNENs]oN*^Att5liwRUArNCOj-l*OUlHXlNCOL«NHO«l,pcSLll^RO* I ;£caHWROAl g 
• OKlMcOLl •rKlNM ti AxIN**OR) .oTRInCUL* tBMNRUNl tSCNACHlHCOLfNRORl 

CONFUTE P 5 EUOU inverse 

A A* 5 * 

00 IS l■•l<UKO« 

AA>ANA(ttl)*«2 ♦ AA 

continue , , 

CC ■ ITOLEtlC*»>l .• ANA^ill . *>A» 

IF 1 ABSiAAl .OT.tCi CO T« 3o ^ . <- ,,-n 

IF first col vector is XErO then first rOr of OUlH* iS 2ER0 
Co 20 IbI.NKOk 
OuTha{1,11»3. 

Continue 

IF First col. is ho-i /Ero, then oo first ru# of at 




10131131 tR 


H 



001 25 

296 

OOlJD 

336 

00] JS 

31* 

ooi Ji 

32* 

ccUJ 

33« 

OCJ3J 

34» 

COlJ't 

3S» 

001 J^ 

J6« 

001 JS 

' 37* 

00136 

38» 

OOIHI 

39* 

0QI92 

NO* 

OQIHH 

4l» 

OOIHS 

42* 

ooisa 

M3* 

031^1 

44* 

COlbH 

4S* 

□ Ol:»S 

46» 

0O1S7 

47« 

Do163 

466 


49* 

Col 61 

50* 

03161 

51* 

03162 

52* 

03163^ 

53* 

O3I 66 

b4» 

C3lo7 

5S« 

03171 

56 » 

0317 2 

57* 

03172 

sa« 

C3172 

59« 

OOl/H 

6j« 

30l7S 

6 1 • 

00176 

62« 

33177 

* 63* 

03230 

64* 

03203 

65* 

03204 

06* 

03206 

67* 

, 032U6 

63* 

(1 , 03ZU6 

69* 

^ 03206 

7J* 

\j 03207 

71* 

^ 00210 

72* 

0321 1 

73* 

03212 

74* 

00215 

75* 

00216 

76* 

00223 

77* 

03221 

7B* 

33 2 22 

79* 

03223 

80* 

03226 

61* 

00227 

82 « 

03231 

63* 

0Q233 

84* 

0023H 

as* 


JO 

HU 

C* 

C<!» 


so 

e»T? 


70 


C* 

C» 


80 


C* 

C* 


as 

c* 

c* 

c* 

VO 


92 

9B 


96 

iOn 


Oo Ht 

fCOLiD - MmUflU 
CONTINUE 
TkaNSPOSE f|«Sl 

thUL T« 4NsP I F Ct. 1 1 U i%L ^ | , nKUW i l l 

miET Tft4t.5POSt «Hif»Jr*AL rOL^t PWCUJOC f IS SCHLAK 

CaXL H0ETHXlTCOt»JCoL,JCAtA«,l,rtROifil»l,<lKoif| 

HOtl SCAtAK BT T^AHSPO^E 

SCAL4R « 

Do 5o I»i*NHOA 
ouTHitii*ij*fcot in • S^AIaB 
Conti NOE 

CONI I HUE 

Oo loo >C«:?»r,COJ 

HaitC* J 
Du 70 

continue 

CAtU KULlHXt 0 UT»tA» A<iDA tN.HHUBi l iHCOL*NhUaI 
Call HOtTHAlAHA »OK»FC^U»HK0S*Hi 1,NNU4 ,hCoi i 

Call s u b h x t a < » f c » ck t * » t »hno^ > 

UETCimlHC ir CK*Z 

AA«J« 

Uo Jo H»i •HR oaf 

AA«CXiM}»*2 ♦ XX, 

Coiitmt 

Cc« i TOLeWc* • 2 I • ANAx 1 1 1 • • X a 1 
iFlAQSlXXl^vsT.rC^ UQ V 3 

If ck » 0 riiEn Co'iPuT*- o< 

Call T«ANsP<oit..iT«, ^,t tsiCou.ii 
c iLL hulTik J ar lA v«i*^fKua» i .iicol^ 

Call HOLTnxiorR ,jk,scai-aK,i ,*i,i , 1 sHCOli 

Oh» 1 , ^ SCALA»t 

uo as a*i . i*to,r 
OKini « bkIiii/dm 
Continue 

To 9S 

IF CK»IIEcn tBX=«C^^ 


CALL TKANSPiCK .TCOL.HaJAt't ttmoJi 1 I 

Call HULtHXlTCaLlCKiSC^LARiI .N«01ttI*l»NRU«) 

SwALar * 1*7SCaL^H 

OQ 92 H»I*|n04 

o<«Hi » Tcolihj • Seal ah 

continue 

CONT JNUt 

Call HOLTHXlOK,tt^>SCHALHiH»1 # NHO* ,r<COL# I ^ 
Call SU8HxlaUrHX»SCRnC«|OuT4X|NfHHOa»«LoL J 

Du 96 H»ltN»<D5ll 
OUTHXIK.KJ * BkI«I 

Continue 

Continue 

KLf URN 

Eho 



C«0 iJF COKPILAtIaM 


NO UlAutUfSTlcS 


e rOK ».HAO¥fl.KAOV 

UmVAC U08 FOKTKAN V Extc II LeVpl. 

This cohhilation »as ooke on ?s har 70 a1 


”IE*EV« level EIEOIOOIO^J 
ISI3>:33 


subroutine lhaov Entry point 5t«.7is 


STORAGE used: CODEtl, 007262{ oAxAlCj L3577‘»J bUaNa COhHOvIZ, Ooo&OO 


external references ibeock, nanei 

C003 PSIKV 
DOO 4 MUL7HX 

0005 SU 6 KX 

0006 TRAHSP 

0007 NNDUS 
ODlO N102S 
OCJI nrous 

00 1 2 NlOIS 

0013 NERR3S 


STDRa&E ASSUNHeNT iBlOCk, type, R»t*T 7 v£ L0C4TION, NaHEI 



OCOQ 

000103 

IGC5F 

nnoo 



ICuIf 



ecqsv? 

IQ^Zf 


QDOD 

ODD 1 1 6 

lUPSF 

P03D 


rroi?** 

10371 

cool 



igJSb 


OOUl 

acTHas 

1 113a 

CGpl 


E-:3S3l 

114Qt> 

3031 


P 3 JS .34 

1 1^3b 


OCOl 

004003 

1Z36G 

coal 



127“l> 

203 1 


C5Q321 

13L 


or»oi 

ocHaaE 

M^ 6 g 

C03I 



1375t 

3001 


o:.ci3| 

lIpG 


CLQl 

0 CCIH 6 

IH^C. 

CtQl 



1 

CCjl 


D546b2 

iHf^b 


OCOi 

CLI6I5 


SPPl 


Z ? C < - * 

1£>3g 

ccoi 


C^5l7l 

iS^CG 


DCOt 

DCI632 

l$Ct 

ocai 



l 6 i 2 W 

2 C 01 


CC5>45^ 

1622 ^, 


DDOl 

0tl6*i7 


cool 


171? 


3001 


CS1722 

173l 


OCOl 

CC2I67 

l«CL 

0 C 20 



1979^ 

0C51 


CS,245l 

200 L 


ODOQ 

. 0CC251 


0003 


1 

22a3F 

02C3t 


C5D32S 

20 pi 


DDOJ 

0CC2S7 


0031 



Zlhb 

3 C 3 1 


DiE4c> 

235g 


000 ] 

OCC****!} 

25pG 

0 C 2 l 



Zbav 

3033 


cciiaj 

250 


0001 

DCG477 

27bQ 

cooo 


2*!C3f.l 

2?V9F 

0001 


ccr32i 

30L 

0000 

0CC371 

y2Z7f 

ncao 


3r3*n^ 


OOPi 


LSC422 

3Q0^f 


0000 

DC0H43 


C0DO 



3r^bF 

0033 


CCC471 

3U0^F 

OOGO 

00 CSI 2 

3al2F 

C030 


zzzsi'j 


00CS 


0 C 2 B 2 ^ 

3Q10F 


CCQO 

CCCSH3 

3G25F 

CDO0 


cri3b6f 

30-^ Of 

5C03 


CitS77 

3a3bF 


GOOD 

CCC617 

3csor 

OCJO 


•2^631 

3^60? 

0023 


c 33633 

JPOPF 


0001 

ttci><>3 

322a 

0031 



326 b 

OCOl 


03*1 127 

330L 


OCOl 

003645 

3t>2(^ 

OCOl 


331661 

362g 

ncel 


C2=7|4 

37 lb 


OOP] 

OCi 106 


OQOf 


3?1 1 J 1 


occi 


CC1IS2 

*fSlb 


COQl 

CC5404 

52Ct 

0031 



521b 

3031 


CL1S7S 

S«flb 


DDOl 

002111 

6«»ca 

0001 


^ 021 1 J 

6b3G 

3C31 


L3t76Z 



CCPt 

GiC 1121 

73L 

GD3I 


33112^ 

7-IL 

CoCl 


Ce25lZ 

^52b 


OCOO R 005000 

k 

cooo 

0 

ZZZZ2^ 

cc 

0DG3 

K 

Cfc3S4t 

CX 


oaoo 1 

I 000204 

HEW 

0030 

R 

tarrsf 


oca3 

H 

CCC242 

HX 


C'OOO 1 

1 DCC{«22 

IIKIC 

0C3C 

1 

3rrc6* 

Ij 

0C23 

1 

D 0 C 221 

IK 


0000 1 

1 -QCODlO 

IttHQ 

8030 

I 

1 ^ 

3 

0Q33 

1 


Jj 


CDOO I OOCOl/ 


0000 

X 

35Cnit 

K8^ 

•*C 0 u 

1 

C5S:o? 

LIHIT 


0000 j 

1 BCCD2S> 

LVZ 

0030 

1 


tP3 

305 

1 

V’jVt2] 

LPH 


0000 1 

1 000032 

tP7 

OOPO 

1 

CCCr33 

LPB 

frou 

1 

0 00365 

« 




O0UQ 

i- 

oooaii 

NHEAO 

ocoo 

I 

30tnS \ 

HrK 

GOOD 

1 

C.CCU6C 

NT* 

0300 

1 

00305* 

NT2 

0000 

i 

OOodS7 

ntJ 

OQOO 

1 

00DD&5 

nth 

0000 

I 

3rt;r 1 2 
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MX 
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OllOO 

!• 

00 too 
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ColOO 

3» 

COIQO 


001 00 

b« 

Cu 1 00 

6* 

CoiQD 

7» 

OolQD 

e* 

ddioo 


OOlOQ 

ID* 

COIDO 

n* 

OUIDO 

! 2» 

GOIOO 

13» 

30101 


coioi 

15* 

CQIOI 

i 6» 

QOlOl 

17» 

Qoro3 

IB* 

00ID5 

iv« 

00107 

2D* 

COl 10 

21* 

Colli 

2 2* 

COl 12 

23* 

OOllH 

ZHm 

oniiH 

25* 

CailH 

26* 

Cni IH 

27* 

OOI IH 

28* 

031 

29* 

00115 

30* 

On 1 1 6 

31* 

Co 1 1 7 

32* 

00117 

33* 

OOl 17 

34* 

□ 01 20 

3b* 

C0122 

36* 

03123 

37- 

0312*4 • 

JB* 

0Q126 

39* 

QOI27 

HQ* 

00l30 

HI* 

03132 

H2* 

C013*4 

H3* 

CO 136 

HH* 

Co 137 

H5« 

C01H2 

H6* 

001 H3 

H7* 

On 1 H6 

HB* 

00IH7 
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r# 

c» 

c* 

c* 

c« 

c* 

c# 

c* 

c* 
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LhAOV is to COkPUlTe AwAltSiS OF ^aHjAhcE FOK j 
UES iGNSrOESlO^S ^ITh DATAi ANU pLSiONb AlTH ^WUTlPtE 

solutions NOH>»aE Eci^A* I 0>^1>, ANAttStS InEluuES ANALTSiS Of 
COVaRIaNcH ANU AhAttSTs oF KEOHEbS I ORt 


inputs - TvEC TO ^OnS} Ak£ ^,0«fv AKKATS AHoSE BASE AUOKESSEs AP£ 

CUHP< T^D in JrAEtlMU 

SToA&GE ALur " 
l> I MF 1 0t4«; ’ 

OUTPUTS* 



Sub ROUT I N r UHA 

• CHV,UV£CtAKH.)( 

• 4i«ST«PARA,uiuV 

• F»OHASTiNP «UW 
Data a/*--~^« 
DhTa LP/»t*/ 
IriTEGER HltH2 
KLAL Krtx 

heal NPT 
Data lih1t/^ 2: 
or hens I OH Tv EC 
•RINV(NW,Nm) iTV 

• ^^TRPt ,NP) tHM 
St uHltX lN?4 rU U . i*. 
•PwRlClNa, HPI iP« 

• 2tRp ( np ) »pP»?Hn 


Hl^ * 0 
te% « 0 

I ZERO * 0 


f.7 I T VeC 1 bVLC iRKX.rrvtCiR I nV.TVEc » feUA • & INV . hIHV i 

TRPfHHA ,CHvEC.RAKR*CHA,CH«Af aTAkS iClHV|tl”X|2«RA, 

. pPb I ^ ZT^P • pPfiRK t LWKK » A AKr f 
M t NU I u alTS » tULEnC • tu»S J 
»/ 

.HtW , TEv 


tUK ^ 
) i G 

J * 

#UP ^ 
P) 


»P.UU> ,avEC(HHJ •H^UlNt1>Uvr3 ,cHVlNw) • 

AiNUtUOJ •GINVIUW.Nm ,HVtC<HNil »WKR«NUaNW) t 
G«<KtHNiNNJ fC^Vtc^NlO *C}VECINH| 

) »HrtX*HP*NU» *LHAlHu,uUJ t4*iKK • 

• H I H V i • N N ) » H I N V TN g • N P J I H N y t N U , W W 1 , 
i»<OHKbTt I ) »«Ahr INWiNR* 


• Read i.t Tvec a^j xmx aT Firsf analysis 
^KITE(6, I07Z) 

I0n2 PoRrtATUHl |///SX» ’ I raPUT DaTAM 
HKLAO«uy/Hl 

I F t H o U C nU i « j , G I . O J H Rt. ■ Hr E aD> 1 
NwKITE * HU -7 . 

U.VRITE » N *filTf/» 

if lMOO«NWHl TE.R J •GT* j)HaK| TE « NftRlTL ♦ i 
iFtNALYS-Grt I )(,0 T L 33 

^RiTEiB. io::3» 

1503 FurtMATl//SAt«Y vEC oK AhU X 3aTKIA») 

Do 2n I«1«.4P 
K«I 

Do 10 Jsi.HKEAn 
• If ♦ 7 

iKlK8*GT*Ng)KB*N^J 




liffliayii 


OOlSi 

ao# 


00157 

ai« 

1000 

0916Q 

52« 

0gi 6 i 


10 

Ooi«3 

58 « 

0g|b6 

5S« 

1001 

Sour 

56* 


. 0qi7D 

. 576 


00172 

sa* 


00201 

5Tt 

2M99 

50202 

606 


aozoH 

6l» 


00205 

62* 


00210 

63« 


0021 1 

68® 


0o2U 

- 6b® 


00221 

66® 


00222 

67® 

12 

00225 

68®. 

1008 

00225 

6V® 

13 

00226 

70® 

20 

00230 

71® 

30 

0023 1 

72® ' 

00233 

73® 

lOOS 

00235 

78® 

00237 

75® 


OOZHO 

76® 


00253 

.77® 


C025<» 

78® 


00256 

79® 


0025 5 

80® 


00255 

61® 

^5 

0C257 

62® 

0C2&2 

83® 


00271 

68® 


00273 

6b® 


00275 

86® 


00277 

67® 


00300 

689 


00302 

69® 


VJS 0p3l0 

90® 


CC311 

91® 

18 

^ 00313 

92® 

IS 

00315 

V3® 

SC 

CC3I6 

98® 

00320 

95® 

icf»r 

00321 

96® 

00325 

97* , . 


00325 

98® 


00330 

99® 


00331 

ICD® 


00333 

101® 


0035 1 

tC2® 


00352 

103® 

60 

0C35 5 

108® 

00357 

lOS® 


00356 

106® 


00360 

107® 



K£ADi5,i0fiai (XHXU |IC(CI 

FoKH4Ti8Fia*iJ 

KmKU 

Coktikue 

Head 15 »i nn i itvec ^ ^ 

Fo««*T(ria«6) 

Kl m J 

ir(HO«LE*Kl }iCl«N^ 

«KlTEf 6,2*|9^j YVECI I j (AhX|1 ,KICI »KK*ltKl } 

FoRhuTI / 7I4X»FIo*6^ I 

lFlHQ*tE*riGO TO O 

Ik » 8 

00 12 J»I*H8H1TE 
IlKK « IK ♦ 7 

Ir i 1 IKK.&T.iiQ) I I^K . h4 

IfH I T€ 1 6 , von H > i I 4 • I kK ) , I K K« I K , 1 1 KK I 

IK • I IKK 

CONTIMUE 

Format t / 1 kx i 8 ( p i D* c , 3x 1 1 

continue 

Continue 

Continue 

•t?ITEl4. IO?5) 

FciHMATI//Sa««t vtC^oR AhD k HaTRIx*! 

Ou SO 

K«1 

Oq no U*l.NHEAn 
K6 • r ♦ 7 

1 r I X 8 . G T • NS? ) 1C B £ N ^ 

Ne XO C 5 » j ocri ( Rmx I 1 » kX I » XK«iC , K8 1 
X»K8 

Continue 

MtAOfS, lOfil ITVfCn I 

AKITEI 6,2S7V J TvECI | ) , C««X i 1 rKKl tKX»J »Kl > 
lFiNS»LE«74tsQ TU IE 

Ik • 8 

Vo M »K^Hlrt 

I IKK » IK ♦ 7 
IMIikk •&T«no)|Ikv ■ 

JKI TE 1 1 » 1 C21 ) { hhX I * , I kK I f I XK»1 K 1 1 IKK I 

Continue 

Continue 

Continue 

JieiTrve».ln?7) 

F0F«AT<//5X, »h yECToR ^ND LANBUa MATKiXt) 

OO 70 I»ljNH 
K*I 

OO 60 U»I,NHEAo 

K8 ■ r ♦ 7 

lFiK8,6T*Na)Ke»NW 

Read 1 5 ff loen v « ghx i v tKKi • kk«k .ko i 

K*K6 

Continue 

RtA0c5*l0r>nHVrc^I I 

• Ri TE f 6.28991 H vEC I 1 1 » C «*nx 1 1 ,KkI • KK» 1 1 K I i 
If lN0«tC«7lG0 it 

Ik • t 


&0361 

106» 


00361 

10?P 


&Q36S 

II a* 


.00367 

111* 


C037S 

112* 


0n376 

I13» 

16 

90100 

UH9 

17 

00101 

Il5« 

^ 70 , 

OClOl 

1 )4» 

C* 

OCIOI 

Il7# 

C* 

geioj 

lie* 

c* 

OCIOI 

1 |V« 

c* 

coioi 

120* 

c* 

03103 

121* 


00101 

I22» 


CC105 

U3X 


CC106 

12H* 


00107 

12S« 


Coho 

I24« 


Ceil 1 

127* 


OollZ 

12BP 


0l>1>3 

I27* 


OOllS 

I3Q* 

71 

Ofliu 

131 • 

25fi9 

00121 

tJ2* 

00121 

133* 


09121 

13H# 


00122 

I35« 


00123 

13 b* 

72 

00121 

137* 


53121 

I3BP 


S012& 

137# 


09125 

mo* 

C* 

00126 

mi* 

09131 

m 2 * 


08131 

m3* 


00136- 

1N4* 


00137 

IBS* 


COHO 

mo* 

73 

. Com 

m7* 

7H 

00112 

1H8» 

75 

Com 

m?* 

76 

Com 

ISD* 

C* 

00116 

151 • 


09116 

152* 

C 

09116 

153* 

c* 

09116 

15*)* 

c 

03112 

155* 


C91&O 

156* 


00153 

157* 


03151 

I5B* 

76 

09156 

157* 


09157 

160* 


09161 

1 oi* 


09162 

162* 


OOI 6 S 

1*3* 


00166 

J6H* 

77 

C0170 

165* 



V \> 14 J • 

IlKlC • IK ^ 7 , 

*R 1 TE 1 6 • I C2/I I ( GH^ U » 1 i » J KK» I » 1 1 Kt^ ) 

IK • I IKK 

CoHTl«U€ 

Continue 

Continue 

Compute n*ghxi ^•IDEN l^Ty Hatnik 
Compute base Ar^oREssts fOn subkoutihe psTnv 

LP1»NH^1 

CP2«NH^I^eP| 

UF3«N0^I^tP2 .. 

CP*t«NH^J^LP3 

CP6»N0^I^tPS 

M*7»WH>1+LP& 

tP8*NH*NC^I 

lrltP8*Lt«Etfm tO '2 
CP7 m 1.1 Mir « CP« 

«KlTC«8t2b27lLp7 

F0KMAT(5X^*»**nlMtHSlU^<S or MATM hoDELS CaceCOS NA&E CKPacI^T Q 

•f CHaoV SUDKOUTlNE •• ANaLYSIS ASKS fOR*t|Sf» H0H£ LUCaTI^^NS THA|# 

•available*) 

hcTurii 

Continue 

Call PSlKv l»<MX.«lNir,i<?1»f<N,TOLENC »KOKr*STcl I , A0RKSTCLP1 ^ i AOKKSTCLfZ) 
^•iiORiCST ILP3) .wOHKSrcLP^I •aUNKST ILPS) • wORKSTI Lf 6 » » AONKSTILP^ I ) 

Call «ULTHKlRi«<v»H|iXtN*^Ri»,0,NHtMQ#UQ»N«) 

00 78 

.Do 75 l•^tNN 
lr(J«r9*lM>0 Tn '3 
NKRlUf ll • - ^KK^UtU 

Co TO 7H , 

•aHU,!) * !• • WKRcUill 

CONTINUE 

Continue 

Continue 

Call HUlTnX <G*4X,^Kn,HHA,NN ,«Q,NaiNNiUO) 

UETERNlNE ir TvEC • 0 AnO hVEC • 0 

KjC - 0* 

00 74 1*1. HN 
KX«AA*HVECf ll••2 
CONTINUE , . 

CC • ITQtENC«*7J .• AMaA|I,,,KX> 

IrlABSlXKl •i»T*CC * <iO T® 7^ 

AX • 0* 

DO 77 I*1,NH 

AX • XX TyCCri**^7 

Continue 

CC * fT0Lr;«CA»7) • aHA^I^vvAX) 




OOH7i 

1 a fc'» 


D0H73 

167* 



I 6d9 


03*! 75 

l64* 

79 

Coh75 

173* 

C* 

03475 

171* 

c* 

004 76 

. 172* 


00477 

173* 


DoSOO 

17H» 


0D;>XJI 

I /5« 


Sq502 

176* 


CG503 

177* 


0S534 

178* 


0q555 

17 V* 


20536 

005X0 

18C* 

181* 


Co5l0 

1B2* 


00513 

163* 

c* 

0051 1 

134* 


305 X Z 

IBS* 


03512 

136* 

c* 

Co5 S Z 

187* “ 

c* 

00512 

188* 

c* 

OG533 

1 8 V* 


C25 I H 

rvc* 


Co 5 i H 

1 9 1 • 

c* 

30515 

iva* 


30516 

1 V3* 


CG5I7 

1V4* 


305 20 

IVS* 


30523 

196* 


OvhZH 

197* 


COb25 

iv^;* 


Oo526 

1 9V* 

so 

3053 3 

2 GO* 


00531 

2Cl* 


3353 2 

2d2» 


00533 

2C3* 


CC533 

2G4* 

c* 

C3S33 

2C5* 

c* 

33535 

d. 0 6 * 


03536 

2D7* 


30537 

2Cd* 


Cos 4 0 

239* 


3054 3 

210* 


00544 

211 * 


CC545 

212* 


C0546 

2 1 3 * 

9C 

CC55C 

21 4* 

00551 

215* 


C0552 

216* 


03553 

21 7* 


00553 

21b* 

c* 

D0555 

2 1 9 • 


0C557 

2Z0* 


CC56I 

221* 


DC56| 

222* 

C* 

0Q561 

223* 

c* 


litHa * I 
<*6 in 18D 


C|ji>^P||T£ •HhX*^HV|iC 

LP3»f|W'^I^tP2 

tF4=«f4Q^l>tP^ 

U*7 

ti 

CAUL^lif^fvVlj^n^lMVjAH^ f 1 1 ,*OHi^br C« PI 3 jWOHnsT iLFi^l 

• j^lOHKSl t tP3J f ILP**3 »* JPK5TILP5»3 tLP«» 3 • i LP^ 1 1 

Call m^lTnxssma 

Call I 


Calculate hhxckjw^ **TyCc“cTvEc 


Call H jlI haIPJ-x ly jtilil 

CALL i I »fm»»i»3l 

Cl»1P*P£ Ch»EC TD «3EC 

^x»o. 

A»Jt*D a 


«x o » 
t?u 1*1 

Xa^CHVEC M ^ • Z ♦ A^ 

«^w««itci 1 1 *♦? ^ 

P A * C ^V w C ^ i • * H * I' C ^ 3 ■ ♦ ♦’ X 

CUUTIAjE 
CX*XX«2»r A* 

HA*XX*Aa 

H|.» TDlEUC»* 7 • »Xl M • tHXl 

if i ASStCXI « 1 


C0HPx?fE Cj^iLC to TVtC 
X A * V • 


f»tf*S« 

ir > • t * 

CO 9r f*3,»H 
Xx"CTVECi I ^**7*K^ 
#/A»TVEC l l J • •2 
?»X»CTvEC 3 1 3*Ty|;C^i3 
CCMXLJE. 


♦ WX 


Cp.^AX-2*^%^ht, 

H*i - TOlEHC-vZ •■ AHaXI t J *,nX) 
If UBStCX3 •LT»kx 3 TLfr * I 


Ir lhro.EC*3 *Am/#i^EC»£y* I 3 gu To ITc 
IF IHEC»E^*3 t Af<0.‘^Ew,tL*C3oO To IS0 
IMHEQ^rEy*r.A^n.Tcy*eW*i3G^ TO ItU 

HEy*;?‘AtSE AliD tC^^FAL^E 


0C563 

00S6B 

00b66 

00B70 

00570 

00570 

00571 

00572 
005V H 
00576 
00576 

00576 

00577 
00600 
00602 
0 D6D<i 
0060M 
006Q*< 
CO605 
0O6O5 
0C6U5 
0C605 
006U6 
00607 
00612 
CC615 
00617 
CC620 
C0621 
00622 
0Q623 
00625 
00627 

00630 

00631 

00632 

00633 
0a&3<< 

00635 

00636 


5o6H1 


OObHl 

3Q&*!2 

SQto43 

S 9643 
06 43 
20643 
30644 
00645 
00646 
30647 
006o2 
Oo6bS 
00656 
0D6S7 
00660 
09662 


224 • 
22B« 
2266 
227« 
2266 
22^* 
2306 
23l« 
232«» 
233* 
234* 
235* 
236* 
2 37* 
23B* 
23 V* 


2H0* 

241* 

24-2* 


243* 
244* 
24B* 
246* 
24 7* 
248* 
249* 
2b0* 
251* 
252* 
253* 
254* 
ZBS** 
2bo» 
257* 
258* 
259* 
260 * 
26i* 
262 * 


2*3* 

2o4* 

26b* 


266 * 

267® 

268 * 


269* 

270* 

271* 

272* 

273* 

274* 

275* 


2/6* 

277* 

2/8* 

279* 

280* 

281* 


>HlTF(6f ZCO-I 
2090 ruKHATl////*5l 
►KITEI6.2003J 


EKr 


*kITEl6.2003j . . , 

2003 F0KH*T«//,5jl,*ftEJtCT M*0> wllH PRObtlrPt.l fHROR 1 «»K05 1 T YPt jl ^ 
•OR>»0.*.//f5»,»NtlTHiR THe ASSUHtO HE5 TrIcTI0N5 NoK THE hTPoTHEsIs 

• IS MATHEHATICalI-T c0n5iStei«t.* » 

, return 

tSo »RlTE(6.2oO;i 

5RlTE16,20Ci| „ . . , 

2C0! F0K«*T(//.5a,*m(C 1 |S * LfOlTlHAlE HTPOiHesIS BUT rHE^MOOLt- IS nOt 

• HATHEMATiCAULV consistent.* i//i5Xi*5Etof‘SlOtR Thf aSSUhEO RESTkIc 

•tions.’i 

return 

160 RKlTE(6.2nG;) 

2002 F0RHATt//5At »Re jtCT Hlgl RlTH PKOblTTPE I ERN0K1 "PrOB ITYP fc 11 tRKo 
,R J»8. » »//.5 A. *thC AsSOHtO RESTRICTIONS A^E HAThttlATlCAUUT CqnSISTe 
6NT. but The hypothesis is not C0HSIS1EMT*»| 

Return 


V 


l7o 


172 

1/3 

174 

175 


IF HMXUMxI tHHXl U*4XHXU*RIMVIR«X^ 1 

ItESsAOE HEtdtf 

CALt HULTMXCRiK,v*HHXt«f»^I<iNw.NMtNQ,NQiNH| 

OU 175 3*1 fJia 
Do 174 I*l»IiQ 
If l3«E0,i JCiO To 172 
• KRU.I) » 2. - nKKcJ.ll 
UU To 173 ^ 

»XRtJ,n » 

CONTINUE 

Continue 

continue 

Call hulT hx i xmx « **kk » hhX • Np • fig • uu • np t 1 

LPi ■ MP ♦ 1 

l^P2 s NP ♦ 1 ♦LPl 

Lp3 ■ no ♦ 1 ♦ L*'^ 

LP5 • NP ♦ I * L^4 

LP6 K 143 * I ♦ L^^ 

LP7 s NF * I ^ 

UP 6 ■ NP ♦ Usi ♦ • * LP^ 

\ju ro 7 1 


• N£* Hrtx 


IFtLPR*ar *L1 11 T I 'jU i u 7 1 ^ . 

Call pSiNV CH^X ,HlNV,i4P»N'l,TOLENC**UKK5T { 1 j ««ORKbT,| Pi) t«ORKST (tP^I 
• fl»ORic5TILP3l *rt\lK^^T :lP*<I • ftOrtKST I LP5 > i*UhKstClP 6) f juOHKi» T C LP7 j | 

Call HULT«X(«IriV»HH 3 iNP»NU,NaiWP J 

Call HUlThA IhHX • itig »NWt rIN tNu » 

COHP 4 RE • IF wul. tQUAL* PRl IT HESbAtiC ARD CONTiMUp 


XX » o« 

44 « 0* 

4X » 9i 

Do 177 I * I .Nu 
UO 176 J * ItflN 

XX ■ AARR(3tl)®»7 ♦ aX 

• HHXtJtn«®2 > 

• X ■ HARfUJtll • HKX1U»|) 

176 CUNTINUE 

177 Continue 


4X 


1 


m 


0966H 
03&6S 
QObhi 
C0667 
50i7| 
0067 J 
036 75 




* 2 a 2 # 

2 tt 3 « 

2:1 il 9 * 

2 d 6 « 

2 a?« 

2 as« 


0067S 

289# 

00876 

29Q« 

D0876 

29 !• 

30876 

292» 

3067 6 

293* , 

006/7 

294» 

3o6 7 7 

2V6* 

006 7 7 

296« 

00700 

29 7# 

C 07 OI 

^ 29d» 

CO/Ol 

299# 

00702 

303* 

C07Q3 

jQi • 

00704 

302« 

007U5 

303* ' 

00 7 06 

304* 

O 07 J 7 

JDS* 

GO? 1 0 

JQ6* 

0q7U 

337* 

007 ! 2 

3oa* 

0071 4 

309* 

007IS 

310* 

007 15 

3J I* 

00717 

3|2* 

00717 

313* 

00717 

314* 

00720 * 

3|S* 

CC72I 

316* 

00722 

3i7* 

C072J 

3i8* 

00/24 

3|9* 

007 27 

32D* 

00730 

321* 

00731 

322* 

C0732 

323* 

00734 

324* 

CS735 

32S* 

Ct?736 

326* 

00737 

^27* 

00741 

328* 

00743 

329* 

OQ7 4S 

33D* 

C0745 

331* 

00745 

332* 

0Q746 

333* 

00746 

'334* 

30746 

335* 

00746 

336* 

00747 

337* 

00747 

33B» 

0o747 

339* 


It9? 

178 

C* 

C» 

C* 

I 80 

ll 


C# 


c« 

c« 


CA • AA « I 2 » ^ 

Ha • XA ^ 4 k; 

M4 ■ irOLrUC*i7)*Afl4XlM*.HAJ 
iFt A8 S(cx> ro i7h 
•RlTEi^iZoiSl 
•HlTEC A, 19791 

FQRHxT</^/bA,«THE H/pOT^E5i 5 rUHCTiUH *UAm«* DA»B£T A i |S HOJ E&TiHABt'E 

• • analysis COf^TlHuHS* • * 

CoNTinue 

Calculate tmen c^hpaRE^hhx ihhaI'*-! • tHvtc»<»HA(H:Hx^9 tvec<*ihyEc-cima® 
(RMX^jTvcr 

Call hultha .hihv .c^^hx ,nn » fiii iNH* H y i 

CHYEC COfiTA I NS ( HVLCwbHX IR I M v ) TVEC I 
Call HULrrlAlCHw A»TYEC.f5'HVec.NH»NH, I 

CALL SUaHXtHVEc tCrtVec.CHVLC.IiNi I f UNI 
Ll>l«NN+l 

UH2«NN^1^lP| 

LP3»N3^1^I P2 
LPH«NN^|4LP3 
LP5*NN41 ♦LHH 
LHd^MQ^i^LPS 
LP7*UN*1^LP6 
LpB«NN*M<I^I ♦LP^ 


ir«LPa»GT.Li«intiu to 7 | 

Call PSlHVlhHX .HUlV.NN»NV|,T0LEHC»A0Rr»Src I I •NOHASrii PI ) tAOHKSrCLP^I 

• »*iOR)csTiLP 3 | .Au«Lbr(LP^I ■aOrKsT « LP^> t AOrPCsx 1 LP 6 > * nOrKST C LP 7 1 ) 

Ir. n Jria 


190 


znoH 


I r LJ ZE*to»Eu • 1 > GO lu 20^ 

Call HOLThaCMHa V»UMv »C*'RA.NN.NW t NN,NN»NyJ 

g^Ra contains 

HvEC COnTaIhS HH'*“*ChvLC 

Call KULTMAU,i,KX»LHVtC»HVrCtNNfNN, l iNNirH) 
AX*»0. 

*A ■ 0 * 

Oo 190 I»I»HW 
Ax»CHVECt 1 1 •*2>xX 
«rt*HvECI 

HX»ChV£C t n •HV pC< H ♦ 

COHT INUfc 
CX*XX-2»Rx^r<ff 

IIX*XX*MT<i 

H* » T0LtNC»»2 •, AM*X1 ‘ l*,HX) 
iFiABSlCXl •LT.v»i^«*0 To 20o 
»RlTE« A,2o?3) 


i^s . 

r-m 


Jar 


WRlTEfA, 20 ^N) . ^ . 

I F0KHAT{//*bX, •RE'^lCt RjTH PRom TYPe | ERKOK > «PR 08 i TYPC jl £Hr 

•OR)mg, t ^//bA* • TKC RESTRICTIONS ANU HtP0T»ESE5 ARE SEPARATELY C0 nS| 
•STENT. BUT NOT jOlHTLY COkP*T IBLE • • J 
RriURM 


i: 

c* 

200 

c» 

e# 


Return 

calculate M.xUtI i*lRHX^|RHX-«HHX4lHHXI 

A»XR CONTAiHSIH*iH 

Call hult« x r« i n v 5 hmx . * ar , , ln>nr , n« . nn i 
I»AK5 contains irRx^trk^ 


|C0750 

34 0* 


|Cq751 

3 H 1 • 


|C07bS 

3H2* 


f 0O7S7 

343t 


f 0076l 

344* 


f 00762 

3H5* 


1 00763 

3H6» 

202 

f C07fcH 

3h7* 

203 

1 0n765 

3HB* 

204 

* Co? 6 7 

3 4 7* 

205 

r 0O771 

350* 

Co772 

3Sl» 


i 00772 

352* 

C* 

i 00772 

353* 

C* 

* D0772 

354* 

c* 

1 00772 

35 b* 

c* 

i 00772 

356* 

c* 

1 Do773 

3 57* 


r 007 7 3 

JbB* 

c* 

C o 7 7 H 

35V* 


^ C0775 

U 60* 


C077b 

361* 

c* 

CC77S 

362* 

c* 

. C0776 

3& 3 * 


. C0776 

364* 

c* 

00776 

365* 

c* 

Q0777 

366* 


0^777 

36?* 

c* 

C0777 

36U* 

c* 

? Do777 

369* 

c* 

' C0777 

370* 

c* 

; 3C777 

37 1* 

c* 

• DlOOD 

37 2* 


; OlOOO 

3 73* 

c* 

' OlUO) 

3 74* 


OlOOi 

33b* 

c* 

01002 

376* 


ntOD2 

377* 

c* 

D1G03 

37 B» 

ClOOH 

377* 


OlOOb 

3B0* 


DlOOb 

3b 1* 

c* 

C l 005 

3B2* 

c* 

0 1 QU5 

3B3» 

c* 

CIQ06 

364* 


0I0D7 



0 1 0 VO 

3tJ6* 


0101 1 

JB7* 


01 0 12 

3UB* 


0IO13 

3H9* 


DlOlH 

390* 


0IQ15 

3V1* 


SISU 

iV2« 

393« 


C1Q2Q 

3V4* 


0102Q 

393* 


01021 

396* 


01022 

397* 



Call I 9 KH> NV i U^ • Hh* 

VC Zfl5 J*l tHQ 
I/O 20H 1*1 •110 
U t j*Eu. 1 1 00 -ru ^Z2 
WKhb{ltJl • 0. - ^KRbll'iJ) 

OC TO 2o3 

WtCKSlltj) « !• - nX(^SU»J) 

CuNT I NUE 

Continue 

Continue . 

C*LL SUbHx < *xf<«; v” rf^tt »N 0 t ua> 

C/VUU MUlTMXIXHX » «P * f^W J 

H*MHx , , 

CAUCULATH C»( I-iKnX+ ) R5'xJ t?TRPC XHX) U« IKmX^ j rHXI 
*»KKS contains 


Take tramsposc amx • Xtrp 
CALL TRaNsP tXMi . A IKpihP .N o.nPsNQI 
WKX CONTAJMS X1K^’*XHX 

Call HULTNX(XTkP»^HX i*f^R»NW,NPiNR,HW»l»Pl 

LPH « NP I ♦ L^-^ 


r!^XH « XTHp I XhX) i I- iRHX'*’ IR hXI 
Call MUlThX t ^K« . ”KHS i E*’HX ,NU , NW t NW>KU tUU i 


C«CHX 

Call mulTux »lukk |Onx ,nw iNw iHW»nQ> 

Calculate: z«yvf c^anx i k^x-*' j t vrc-xrtAiHp»x^ i ihvec-ohx tvE.^i 
CHVEC contains (HVtc-0*^XlRRx»)TVEC) 

CHV COUTAlNStHHX'^ICrivE.C 

CAUL HULTHXtHM.V*CHvtC»CNv,NUtNN» I *NRtNN* 

ZaHKk^XMX • CHV 

Call MULTMXC xHX.ChV .Xu^mNP »n«t l 9HP*UUI 
CHV NO^ CuUrAlfjS ^RhX^/*IvEC • 

Call MULlHXlRINVf l VtC|CHV,NQ,tlrt. I tHttiNH) 

ittST • XMx*CHV. 

Call HULTHXtXHXtC»lV,Xt4^1 *NP.NUf 1 iNPtNwI 
Call SUbHx^YVLc •^’»bTi ^“ST*NP. l iHPI 
Call SUHHxCXaVI ••-»«HKiX’’KX,r4p, I ,NP| 

2«XArK 


calculate riP-Ki - lRn*xnxtCHX^)XTKP • SOlVE FoH ulIRt^UNOEO ihT*) 

LPl« NU ^ I 

LP2« NQ ♦ I ♦ iPl 

LP3 » NO ♦ I *•- L*^^ 

LPH * NO ♦ I ♦ L^‘3 
LP5 * NQ + I ♦ 

LFtr B NO ♦ i L^^ 

LP7 B NO <«> I ♦ L^*^ 

LPB B NO •NO ♦ I ♦ Lp7 

iT<LPA*GT*Llni r>V^ 70 

Call PSlUv(CHX.ciNV ,HU»NW,T0LtNC>AUniC5l I I J .^OKkOHi pi I 

• iWORkST {Lp2| .WOK<bT(LP-^ l » aUOKSI I LPH ) t •fURKbl I LPb ) t AURK^T I LP^» ) 

• *00RKST(Lp/n ^ 

Call hul^hxicimv »x fopg^^ jf.RK«MOfWWpNpfNutN'<) 

CALL MUlTmXIXHX ,Fwhx iP^ST ,NPiUotNP »NPtNuJ 




iMMHUMai 


0 


022 

3VB« 

0 


022 

39V • 

0 


022 

H00« 

0 


023 

HOl» 

c 


029 

HD2« 

0 


027 

i03A 

0 


030 

HOH* 

0 


032 

HQ5» 

0 


032 

HQa* 

0 


032 

H07« 

0 


032 

Hob* 

0 


0 33 

**ov* 

0 


03H 

HIO* 

01 

035 

Hi i* 

0 


036 

H|2* 

8 


037 

Hi3* 

0 


OHO 

HlH* 

0 


OHI 

HI3* 

0 


092 

HI 6* 

0 


093 

HI 7* 

3 


0H5 

Hid* 

c 


OHS 

HI 9* 

0 


□ H6 

H20* 

8 


OH7 

H^2l* 

8 


050 

H2Z* 

8 


051 

H23* 

8 


OSH 

H2H* 

8 


055 

H25* 

0 


0P7 

H26* 

8 


0S7 

H27* 

9 


05 7 

H2d» 

0 


0:>7 

H2V* 

8 


OoO 

H30* 

0 


06 1 

H31* 

8 


062 

H32* 

9 


062 

H33* 

8 


06 2 

H3H* 

0 


062 

H3S* 

0 


063 

H36* 

8 


063 

H37* 

8 


OoH 

43 d* 

0 


Q65 

H3V* 

0 


06$ 

hho* 

c 


066 

HHl* 

0 


070 

HH2* 

a 


Q7l 

HH3* 

8 


0^3 

HHH* 

0 


105 

hhSa 

0 


102 

HH6* 

0 


103 

Hh7* 

3 


105 

HHd* 

3 


105 

H49* 

3 


106 

H50* 

0 


no 

Hbl* 

0 


1 1 1 

H52* 

3 


117 

H53* 

0 


123 

H5H* 

0 


I2H 

H55* 


f 

c* 

c* 


2IS 

€• 

C» 

C* 


22ft 

C* 

c* 

c* 


c 

c* 

c® 

c 

c* 

299? 

3301 

3005 

3S02 

3003 

3009 


P»5T COHTAr«S lJt”AlclNTMTRl»J 

subtract IULkTiT^ HATH^X rKO« PAST AfID TRACk 

TRTOT-0. 

UO 2|S l«lfNP 

rNtTOT • TRTor ♦ U. - PtsTCl.lll 

COMTlHUk 

Ki»FtOATtsPJ- TrTUT ♦ *5 

Compute Rr-«2*TR4AMiiciwV)xrMP-HMxcMU4tf a • souve Fuq r3(roumo ihti 


tPI«RP>l 

tP2*MP-M^LPl 

U»3*Rq^|^LP2 

LPS»WP^I ♦cP9 
tPA*U0^|^LPS 
LP7«HP>l'fLP6 

LPB*NP»N«>| ♦ I p7 
IfUpR.GT.LIMtT) ,<iO TO 7l 

Call ps*‘**’' * 

••AORkST 


>I«GF*LIMtTI .OO TO 71 

PS|NVl4HX.Hl«V,NP»Nt4,TOLEWCt*»ORKSn l|,iiOKitST|, PI J .hOKKSt 
ST ( LP 3 } • A DRXiiT I cP*! 1 • AUK ST i tPS ) t rtORKSf I LPA I * J C LF ^ I I 

MUt Tma CmH> .rtlMV.ppMK.^.iiP.Mij .Np .IIP.mWi 


0KKstCLP2| 


Call HULTH*(HHJi,rtIH¥.P'^«f<iC,«P,W«4,HP,llP,N«l 
Call su8H*«PASt .*’P aba.,ppahii,np»»ip.mp> 

T r n -V • n 


TltrOT»Oc 
00 270 I»|.%p 
TMToT*TRToT ♦ PP^KKil.ll 

Continue 

H2 ■ FtOATtNl) - THtOT ♦ 


Calculate ot^zt^pixha iLriX^i 
PmST COmTaX is XHAlCuA^ixTRp 
Call HUcTfiX I P AST • XfIliN , X^St .I|P ,HP • I tNP tHp 1 

Call TKANsPi2AKK.»xrRP,»*p, j ,np,ii 
Call mulTmxiztkp»xnsi .' rf, i ,up. ; , i ,Npi 

Calculate i j zmx 

PPARK CONtAI^IS HrtXtHUXjl 

Call MatTHXinMA,*UNVfPFRRK,NP,N4.NP*UP»MNj 


Call HULTaXIPP i^AHK^XASl .MPiNri 1 fliPiuFi 
Call MULTuX l ZTuP »Xh*^ST l » I inP • t i l •«? ! 

Print aov tarlf 

R»XTEIo*27V9| 

FoRHATS iHi • AtiV lAaLE*) 

PkITEI 6,307:i 
MHlTei6*3ori ) I Ii ,u»l^2l ) 

» oRmaT (/» ix,?t 

TOKMATC////! 

IIRITFIP.335Z) 


RRlTFlP,3a5Z) . V „ 

IOKMxTIT6,»bOUPCt PF» .*2 ?, •degree 0F*tT5?,«SUM OF ^RUARES * • Tfi3 i 
•mean SQUARE**ri*»T.*>Ai-UE UF F-STaT I ST I C * I 

AhI TE(Ot3t}?J| 

I0RHATIT6, •VARl A^ lUy* .TaO, t|:H£eOUM* I 
»Ri TEH* •30511 CA. • 


WK1TEI6.3079IRI ,0 
F0RHATt///i1 A.«oUE 
*HlTEl4.305bVR2.N| 


ro 


Bf 


• T33« 13 tTbb,E|S»9| 


Dll 30 

456a 

3005 

Ol I3t 

457a 


Ql 133 

4btiA 

3006 

01133 

4bVA 

C* 

Ol 133 

4 60 A 

Ca 

Oi |33 

4 6 1 A 

c* 

Ql 133 

.462a 

Ca 

01133 

46 3 A 

Ca 

Cl 133 

464 a 

C* 

0113H 

465a 


01135 

466A 


C t |3A 

4o7a 


01137 

465A 


0 H H 2 

469a 


Cl 1H5 

M7QA 


Ol IH7 

.471 A 


Or ibO 

4 72a 


0 1 1 b 1 

47 3 A^ 

250 

01152 

47 4 a 

251 

01153 

4 7b# 

252 

C l I bS 

476A 

2b3 

01157 

477a * 


OllTip 

47tiA 


01161 

47Va 


0116 2 

HBO* 


Cl 163 

4BIA 


C 1 1 6 4 

4B2» 


C 1 i 65 

Sb3A 


01 166 

4a4 • 


Ol 17 4 

4B5a 

3C07 

01174 

4b6» 


21175 

4B7A 


01177 

46BA 

3CI 1 

012G0 

4b9A 


01205 

4VO* 

3008 

0I2D5 

4VI A 

CA 

01205 

49 2 a 

CA 

Ol 2U6 

493a 


Cl 207 

494 a 


01215 

4 95 A 


Cl 221 

4 96 A 

30P9 

31222 

49 7a 


01222 

496a 

C* 

01222 

499 • 

c* 

01230 

5D0a 


01230 

501 A 


01231 

bo2A 


0 1232 

b03A 


21234 

504a 

3010 

Cl 234 

5l5a 


01235 

b06A 


O f 240 

b07A 


01243 

bD6A 

3p52 

01244 

bOVA 

260 

01244 

blOA 


Cl244 

5n A 

C 

01246 

bl2A 


01247 

513# 



f Of<HATr///it6,«l>^t TO ®(N)»iT33i li 

%KiTE(6•^0^&^ 

f^OHHAT CT6 , • i UNAO'^» ) • J 


COhPUTE 


«3»ZTK>*U-XhAiCmx^)xTRPIZHA AhD 
Wl«/lR>'<X«Xl^rtX^|jtTRP - HMX«MhX^IZRX 


X«XlCHX>l tTRP PR^t aNU HHXIhHX^I %S ih PPHiftK 

Calc suaHxiPrtsi •PPnRXi^P^RK.KPiNP.KP) 

Call «UtTMXlPPivRX.,/wKlC »ZAsT*fiP tRP* I tNPft|P| 

CALt MULTHX(21wPf ZAsr •«! » I »«P»l » I •WPJ 
uo 253 I»»l .up 
Oo 252 

IP I J.E‘4« I ] CaU Tn ^^>5 
PrSTiJtl) « 0* - PRsTi^ill 
uo Tn 25 t . , 

PASTtJ.n * I* - n.ST|Jil> 

CCMTIMUE 

C0^*TIN0E 

CliNTiNUE 

Caul MULTHXfiF%<;T*2ifRk t^A5T»»jP*«Pt I »NPtNpl 

CALL MUlThX (ZTKP»C>liSTt'<i I ,NP# I 1 1 iNPl 

Nr v« hp - Rl 

NKH«ri*R2 

F*tt I /Q • mR/kKr 

S|»3|/NRR 

S2»w/Nh 

ARI TEl&.30::7|NKRt«2 ,S1 'F 

P0>iH*Tl/7V.l6.»0^t TO 13 ' • T5^t£ iS. V • 1 h li # E 15* ^ • 1 1 

• El5.Pi 

»FilT£i6.3nin 

F ORUaTI T6, • 1 *DJ) • ) 

5»«nEl6i30TtilhH*.W»V7 ^ 

Format l ///• T5 . •eRKUR* *^ 33 * 1 3 *1 S6»LjB.V^ jei ,ei5.pj 

Calculate u»zTKP«i.MAi 

CALL MULThX I ZTkP »AilRA,'<, 1 .KPt I 1 I ,UPI 
► Ri TE< 6,3rr 11 I ,21 * 

V.RlTEi6,30?9|f^R,»;|f 

F ORHaT I TA , • tot aL •T34» *3 »t 56 tEi &• V ) 

6KlTEC6,3n2i) t2l^ 

PfiKK cOUTAlUSt ) 

Call psinv \ amx ,pokk,,*4p •u(; ,tole^’‘C i50Rkbtc 1 1 ,*oraS 1 (tpi i •aokkstilp2j 

• *«1;RaST (LP ,|5uRK5I (lP**! i aoRiCST *LPb) > ^<UkX 3T 1LP6I twOKK^T (LP' 1 I 
CfcLL MOUTflXiPARAtTVEC.^^VEr .NO.NPi I .NOiHpl 

Form Alt |H 1^1 ///bX FHL ^ESt LImEAR ESTIHa^E IBLE» UMKESTHlCltD is*» 

• •/Z) 

L'o 2 A 0 l*lrlTO 

WHlTE«6,3al21‘>VECU I 

FuKMATIT 25 .tlb. 9 l 

Continue 

Calculate covaaiLncc matr^^ 

Call mllTm x i x trf » am a * w ar • nw . np i no . no • np i 

LP1» MM ♦ I 


j isu in 

■ 0 . - POST Hi l l 

« I. - Pi.ST|J|lj 


Q 


OI 2 SO 



0I2SI 

&1S9 


0i2B2 
O l 2S3 
3 i 2SR 
Ol2bS 
0I2S6 

Bl6« 

i»l7* 

ill 

Si ?• 


0I2S7 

1 • 


0I2S7 

S22« 


01240 

&23* 


0l2bl 

52^* 

330 

01242 

U2S« 

01263 

S26« 


0«24<« 



01246 

* 523» 

3H0 

-01267 


fll270 



01272 

531* 


Ol2/*l 

b32» 

301b 

0127S 

bj3» 

0|3U<< 

b3S* 

3016 

SI30S 

b3b» 

01304 

b36* 


0I3U7 

bj7« 


01312 

b3a« 


01316 

b39« 


01320 

bHU* 


01322 

b4|« 


01332 

bH2» 

30?0 

0] 333 

bH3* 

3021 

0l33<i 

b*lH# 

01334 

bHb» 

3S9 

0I3S0 

bS6» 

0I3H2 

bH7# 


0I3N2 

bMH* 

c 

01312 

SS9* 

(• 

01311 

AbbO* 


01315 

bbl* 


Ol 350 

bS2* 


01351 

bb3# 

3Sl 

0l353 

bbS* 

0 1 351 

bbS* 


01351 

SS6* 

c* 

01351 

b57» 

c# 

01355 

bbSn 


01356 

bb?« 


Cl 357 

b6Q» 


01357 

b6l» 

c# 

SI3O0 

b62* 


-01361 

b63« 


01362 

b6M# 


01363 

bob# 


01361 

bbb# 

“ 

01367 

b67# 

302S 

01367 

b6B» 

1 

01370 

b#V# 


01371 

b/0» 


01372 

b7M 



tpZm ♦ I ♦ LP> 

LP3 ■ N 8 «- t 

UPH .NO ♦ I ♦ lE J 

LI'S ■ NO «■ 1 ♦ cP** 

LP6 . Ntt * t ♦ lP& 

LP7 . NO « I « lP* 

urtt ■ no«no ♦ I ♦ lp7 

Call PSINVtnKH.4i{(.KStN0*N0,T0Lt;NC(*URN5T| 1 > , KURlST 1 1 PI ) iWURKST C tP2) 
««hOKkST(LP 3) ,*0KKST(LP^I ittURKSTlLPSI t«ORK5TtLP«i whORKSTILP^I I 
NTS-0 
MT2.0 
tir3.NT>i*^l 
HT'I.NT3*ly 
iP(NT<l<bT.MtJINTH*rui 
NT l-NT2*l 
HT2-NT1+5 

IF iMT2t&T.NwlNT2*NQ 
WBITEI4.30IS> 

F0RNATJ///7.S* . ‘OLE OMReStR 1 CTES) COVARIeHcE HAIRIAiI 
■KlTCf 4>3t}|tf I (i.pj IJ, IJ*nT I ,nT2» 

FoHHATI//,TI?.a ‘aI •I*.*I*.I5AII 

II • NTI 

JJ«NTI-J 

DO 3SO l•NT^,NT*«. . . 

lFIII.NE*ll4RUEU»3a2llI 

Ir 1 1 I .EO* I I Jj a JJ * I 
IF I JJ,CiT»NT2l JJ»'«r2 _ . 

IFII t.EO*l I-.RI rE(t-*3J2S| I .C-lKRSdl iJItJ.Nri.JJI 

FORmATI/2k.* I • . l2.* I »,* l3A»El!>*9il 
F0R«ATC/2At • C • . l2,» » • I 
lriii.EC>iii|.ii ♦ I 

Cdhtinuc . ' 

IFInT 2 .LT,N<| 6 o lo 313 
Ir <NTR*LT*ri 4 ) 6 r- To 333 

Calculate mean qP tvec 

ThEan • C. , 

Do 3SI N*|#«P 

FNEAN ■ TnEmII ♦ TKEcm* 

continue 

THtAN “ TNE«U / Il«AH”P* 

HPT ■ FLUaTINP) • lf!1tAN**2» 

Calculate Sjoawe*^ H(j«.t*pLf cURREtATlorf c**efficie«t 

Call HULTrtAHTKP*TVtc*CriV,i 40 ,HP,| ,Ha»NP) 

Call hulThaihveclhviMm »nj.i •! »noi 

R ■ R - Hpy 

Call thahspitfec'^iisi fiPM • nPiI) 

Call nultnai2wst*tvlc»02> I »NP>i it »wP* 

42 ■ 02 • net 
R • R / 42 
RRlTEI6,332blR 

FOKHATHH////,lXs«^bLE 4nHcSTRICTE0 SOOArEd MOLfIPLr CORHELAtION Lo 
.EFFI cIENTt ,,//.T3j.tlSrVl 
Call NULTHAlRlrtV^TVcCtVHV .NOtNNi I •NUiNH) 
call HULTHAlAHX.CHV,2«HKtNP.N0t I iNP.no I 
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IF ill .EO>l i JJ « JJ ♦ I 
If t JJ.&T«»lT2) JJ*NT2 

I M 1 1 . EO • I > *R 1 T E ‘ 6 • 3020 » I , r c I N V 1 1 1 1 J» t J-Hti , J jI 
I F « II»EW»U 1 l*t I ♦ I 

Continue 

If «NT?.LT,« wjGn To S 20 
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